Accepted in ApJ 

Preprint typeset using I4TgX style emulateapj v. 1 1/26/03 



GRAVITATIONAL LENSING MAGNIFICATION WITHOUT MULTIPLE IMAGING 

Charles R. Keeton 1,2,3 , Michael Kuhlen 4 , & Zoltan Haiman 5 

Accepted in ApJ 

ABSTRACT 

Gravitational lensing can amplify the apparent brightness of distant sources. Images that are highly magni- 
fied are often part of multiply-imaged systems, but we consider the possibility of having large magnifications 
without additional detectable images. In rare but non-negligible situations, lensing can produce a singly highly 
magnified image; this phenomenon is mainly associated with massive cluster-scale halos (> 10 13 5 M Q ). Al- 
ternatively, lensing can produce multiply-imaged systems in which the extra images are either unresolved or 
too faint to be detectable. This phenomenon is dominated by galaxies and lower-mass halos (< 10 12 M Q ), and 
is very sensitive to the inner density profile of the halos. Although we study the general problem, we cus- 
tomize our calculatio ns to four quasars at redshift z ~ 6 in the Sloan Digital Sky Survey (SDSS), for which 
iRichards et alJ (|2004) have ruled out the presence of extra images down to an image splitting of A9 = 0."3 and 
a flux ratio of / = 0.01. We predict that 9-29% of all z ~ 6 quasars that are magnified by a factor of fj, > 10 
would lack detectable extra images, with 5-10% being true singly-imaged systems. The maximum of 29% is 
reached only in the unlikely event that all low-mass (< 10 10 M Q ) halos have highly concentrated (isothermal) 
profiles. In more realistic models where dwarf halos have flatter (NFW) inner profiles, the maximum probabil- 
ity is ~ 10%. We conclude that the probability that all four SDSS quasars are magnified by a factor of 10 is 
< 10~ 4 . The only escape from this conclusion is if there are many (>10) multiply-imaged z ~ 6 quasars in the 
SDSS database that have not yet been identified, which seems unlikely. In other words, lensing cannot explain 
the brightnesses of the z ~ 6 quasars, and models that invoke lensing to avoid having billion-M black holes in 
the young universe are not viable. 

Subject headings: gravitational lensing — cosmology: theory — quasars: general 



1. INTRODUCTION 

Two conspicuous aspects of gravitational lensing are the 
ability to (i) produce multiple images, and (ii) magnify the ap- 
parent brightness of a distant background source. Large mag- 
nifications generally require that the projected density along 
the line of sight be of order the critical surface density for 
lensing, which in turn implies a precise alignment of the ob- 
server, the lens, and the source. The same condition generally 
leads to the production of multiple images, so in most cases 
one expects highly magnified objects to have at least one com- 
panion lensed image. 

The connection between magnification and multiple imag- 
ing can be important in a variety of contexts. An important 
example, which serves as the main motivation for this paper, 
is the recent discovery of bright quasars at redshifts as high as 
z ~ 6 in the Sloan Digital Sky Survey (SDSS; see lFan etall 
120001 120 ll 12003 ft. If thes e quasars are not lensed or beamed 
(see lHaiman & Cenl2002l and lWiUott. McLure. & ,Tarvisl2003l 
for arguments to justify both assumptions), they are inferred 
to be very luminous (Mb ~ -27). Assuming further that they 
shine at the Eddington limit of their resident black holes (BH), 
these BHs must have masses of several x 10 9 M Q . 

Having such massive BHs at such an early stage in the evo- 
lution of the universe presents a challenge to models where 
massive BHs grow mainly by gas accretion that is itself Ed- 
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dington limited ( lHaiman & LoeblfeOOll) . Even in the context 
of hierarchical structure formation models, where mergers of 
several BHs can contribute to the build-up of the mass, the 
initial seeds are req uired to be present as early as z > 15 
jHaiman & Loeb 2001). The problem is exacerbated if BH- 
BH mergers result in the ejection of BHs from the shallow 
dark matter potential wells at high redshift because of a large 
recoil following the emission of gravitational waves . In a 
recent model that includes this effect, Haiman (2004) found 
that BHs can grow by mergers and accretion to at most a few 
x 10 8 Mq by redshift z = 6.4 without a super-Eddington phase 
— a short-fall by a factor of ~ 10 relative to the BH masses 
inferred from observations. 

If the high-redshift quasars were magnified by gravitational 
lensing by a factor of /i> 10, this could alleviate the need 
for a super-Eddington growth phase to explain such massive 
early BHs (since the inferred mass scales as iT l under the 
assumption that the quasar is shining at the Eddington lumi- 
nosity). Although the lensing optical depth along a random 
line of sight to z ~ 6 is known to be small (~ 10~ 3 ; e.g., 
Kochanek 1998; Barkana & Loeb 2000), magnification bias 
can significantly boost the probability of strong lensing in a 
real, flux-limited survey. If the intrinsic (unlensed) quasar lu- 
minosity function at z ~ 6 is steep, and/or it extends to faint 
magnitudes, the probability of strong lensin g for the SDSS 
quasars could even be of order unity ( Comerford et alJl2002t 
iWvifhe & Loebll2"002l) . However, for a population of isother- 
mal sphere lenses, all magnifications /i > 2 are associated 
with multiple imaging, and in most cases the angular separa- 
tion b etween the images is more than 0."3 ( Co merford et alJ 
120021) . Recent HST observations of the highest redshift 
quasars have shown no evidence for additional images of 
any of the z ~ 6 sou rces down to a splitting angle of 0."3 
iRichards et all2 004). which effectively rules out the hypoth- 
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esis that the quasars are all highly magnified by isothermal 
sphere lenses. 

The obvious question is whether the SDSS quasars could 
be magnified by lenses with a more complicated (and in- 
deed more realistic) lens potential , without pro ducing mul- 
tiple detectable images. Wvithe & Loeb (2002) showed that 
microlensing by stars within lens galaxies can permit mag- 
nifications as high as fj, ~ 10 for singly-imaged quasars, but 
the probability for /i > 2 is still very low (<0.5% even when 
magnification bias is included). Another question is whether 
departures from spherical symmetry significantly affect the 
results. One goal of this paper is to study ellipticity in the lens 
galaxy and tidal shear from objects along the line of sight, 
both of whi ch are common in observed multiply-imaged sys- 
tems (e.g.. iKeeton. Kochanek. & Seliakl 1 19971 IWitt & Mao 
119971 lHolder& Schechterl 120031) . Ellipficitv and shear are 
known to modify the full magnification distribution for 
multiply-imaged sources (e.g., Blandfor d & Kochaneldll987t 
iFinch et all 120021 iHuterer et alJ 120041) : but what happens 
when we restrict attention to sources without multiple de- 
tectable images is not known. More generally, our goal is 
to present a thorough and general study of the magnifications 
that can be produced by lenses with different radial profiles 
and angular shapes, without creating multiple detectable im- 
ages. We simultaneously consider both true singly-imaged 
systems, as well as multiply-imaged configurations where the 
images are too close to be resolved or the extra images are 
below some reasonable detection threshold. 

In this thorough but technical study, let us not lose sight of 
the bottom line: We find that the maximum probability that 
a z ~ 6 quasar is magnified by at least a factor of 10 without 
having a second detectable image is 29% (see Table0. More- 
over, this maximum is reached only in the unlikely event that 
halos down to arbitrarily low mass have highly concentrated 
(singular isothermal sphere) profiles. In more realistic mod- 
els of the lensing population, where dwarf halos have flatter 
(NFW) inner profiles, the maximum probability is ~ 10%. 

While our analysis is specifically prompted by the SDSS 
quasars, it should be applicable to other objects for which 
significant lensing amplification would be important, such as 
high-redshift galaxies discovered in "blank" fields 6 (i.e., fields 
not specifica lly chosen for the presence of a massive clus- 
ter lens;e.g. iRhoads et alJ l20od iRhoads & Malhotra| | 2001[ 
Steidel et alJ 120031 iBouwens et alJ 120031 |Ouchi etalJl2C 
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This paper is organized as follows. In § 2, we review 
the relevant lens theory and summarize our calculation meth- 
ods. In the technical core of the paper (§§ 3-4), we study 
the magnification properties of simple but useful lens poten- 
tials: ellipsoidal isothermal and NFW halos, with external 
tidal shear. The idea is to identify general features and un- 
derstand the parameter dependences. Building on this foun- 
dation, in § 5 we compute the probability of magnification 
without multiple imaging for a realistic lens population, and 
discuss the implications for the SDSS quasars. Finally, in 
§ 6 we summarize our results. Throughout the paper, we as- 
sume a ACDM cosmology with VLm = 0.3, ft\ = 0.7, a% = 0.9, 
and Ho = 70 km/s/Mpc, co nsistent with the recent results from 
WMAP dSnergel et alJ2003l) . 

6 Magnification of high-redshift sourc es by foreground cluster lenses is 
certainly interesting and important (e.g., Hu et al. 2002; Kneib et al. 2004; 



IPello et al. 2004). However, deliberate selection of cluster fields makes the 
probability analysis completely different from what we study here. 



2. COMPUTATION METHODS 

The lensing properties of a system can be derived from the 
lens potential </>, which is given by the solution to the 2-d 
Poisson equation V 2 </> = 2k. Here k = E/S cr i t is the sur- 
face mass density in units of the critical density for lens- 
ing, E crit = (c 2 D 0S )/(4irGD iDi s ), where D ol , D os , and D ls are 
angular diam eter distances between the obser ver, lens, and 
source. (See S chneider. Ehlers. & Falcolll992l for a full dis- 
cussion of lens theory.) The relation between the position x 
of an image and the position u of the corresponding source is 
given by the lens equation, 



u = x-V(/)(x). 
The magnification of an image at position x is 
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The lensing critical curves are curves in the image plane 
where the magnification is formally infinite, and the caustics 
are the corresponding curves in the source plane. Sources that 
lie outside the caustics are singly imaged, while sources in- 
side the ca ustics have mu ltiple images. We use the gravlens 
software b v lKeetoa d200 1 al) to find the caustics, solve the lens 
equation, and compute the image magnifications for the vari- 
ous models we consider. 

We are interested in the cross section for having a magni- 
fication (of a single image, or a combination of unresolved 
images) larger than /i, which we compute with Monte Carlo 
simulations. For singly-imaged sources, we set a minimum 
magnification /i m ; n of interest, typically /i m ; n = 1.5, and find 
the smallest circle in the image plane such that all images 
outside the circle have /i < /i m i n ; this ensures that all of the 
images of interest lie inside the circle. We then pick ~ 10 6 
random image positions in the circle. The cross section for 
producing a singly-imaged system with magnification greater 
than fi can be written as 



' S in g (M)= / du= [ — !— 



dx . 



(3) 



The first integral is over all source positions u where there is 
only one image and it has magnification greater than /i. The 
second integral is over the corresponding image positions, and 
the equality holds because pT l = \du/dx\ is the Jacobian of 
the transformation between the image and source planes. In 
other words, the cumulative singly-imaged magnification dis- 
tribution can be computed by simply summing the images, 
weighted by their inverse magnifications. 

For multiply-imaged sources, we find the smallest circle en- 
closing the caustics, then pick ~ 10 6 random sources in this 
circle and solve the lens equation to find the image configu- 
rations. Now we throw away the singly-imaged systems, and 
use the multiply-imaged systems to compute the cross sec- 
tion for having a magnification greater than /j, but no extra 
detectable images, 



A mu it(At)= / S(u)du. 

'm(u)>M 



(4) 



This integral spans the multiply-imaged region, but the func- 
tion S(u) selects source positions that produce lenses with 
only a single detectable image; specifically, S(u) is 1 if the 
additional images are undetectable (either because they are 
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too faint, or too close to the brightest image), and other- 
wise. In general, S(u) will depend on the specific source, in- 
strument, and observational conditions. Hereafter we refer 
to it as the "single-detectable-image criterion" (SDIC). In re- 
cent H ST images of the z « 6 SDSS quasars, Richar d" et al.l 
( 2004) were able to rule out the presence of extra images with 
a flux ratio relative to the quasar of / > 0.01 down to a sepa- 
ration A9 > 0."3, or brighter than / > 0. 1 down to A9 > 0." 1 . 
We consider both of these SDICs in our analysis. In config- 
urations where there are multiple images that would not be 
resolved, we include all of them in the net magnification. 

3. ISOTHERMAL HALOS 

The isothermal ellipsoid is a simple but surprisingly useful 
model for studying lensing by galaxies. In this section we 
delineate the situations in which isothermal halos can produce 
magnification without detectable multiple images. 

3.1. Definitions of Ellipticity and Shear 

Early-type galaxies, which dominate the lensing opti- 
cal depth at image separations A6 < 4", appear to have 
nearly-isothermal profiles based on evidence from strong and 
weak lensing, stel lar dynamics, sa t ellite kinematics, and X- 
ray studi e s (e.g ., jFabbianolll989t IZaritskv & White! 1 19941 



Rixetal.1 119971 iGerhard et all 1200 It iMcKav et alJ 120021 
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2003b: ISheldon et al.l2004l) . The 3-d density p <x r~ 2 and pro- 
jected surface mass density k oc R~ l correspond to a flat rota- 
tion curve or velocity dispersion profile, and a deflection angle 
that is independent of impact parameter. 

The lensing pro perties of a singular is othermal sphere are 
very simple (e.g., Jgchneider et al. 1992). The surface mass 
density is n = R e [ n /(2R) where R e \ n is the Einstein radius, and 
the lensing potential is <\> = R e i n R- A source at radius u > 7? e in 
behind the lens produces a single image at radius R = R e i n + u 
which has magnification [i = R/(R-R e i n ) = l+R e i n /u. A 
source at radius u < /? e in produces two images at radii R± = 
R & i n i u on opposite sides of the lens galaxy, which have mag- 
nifications /i± =R± /(R±-R em ) = 1 iRem/u (where a negative 
magnification means that the image is parity reversed). 

An isothermal ellipsoid has a projected surface mass den- 
sity of 



l+q 2 



(l + 9 2 )-(l-r)cos26> 



1/2 



(5) 



where q < 1 is the axis ratio, so the ellipticity is e = 1 — q, 
and (R,0) are polar coordinates centered on the lens galaxy. 
(Without loss of generality, we are working in coordinates 
aligned with the major axis of the galaxy.) The lensing proper- 
ties of an isothermal ellipsoid are given bv | Ka ssiolg_& Kovnerl 
(1993), Kormann. Schneider. & Bartelmann (1994), and 
Keeton & Kochanek ( 1998). For a spherical galaxy the pa- 
rameter b equals the Einstein radiu s, while for a nonsp herical 
galaxy we can relate them by (see lHuterer et all 2004) 



R, 



[2(l+q- 2 )] ]/2 K(l+q- 2 ) 



(6) 



where K(x) is the elliptic integral of the first kind. 

Gravitational tidal shear, produced by objects near the 
main lens halo or projected along the line of sight, can 
increase the probability for high magnifications. Shear 
is expected to be common, based on both analytic esti- 
mates and numerical simulations (e.g.. iKeeton et all Il997t 




arcseconds 

FIG. 1 . — Source plane for an isothermal ellipsoid with ellipticity e = 0.5 
and Einstein radius R em = 1"0. The curves show the caustics. The small 
(large) points outside the caustics indicate singly-imaged sources with mag- 
nification fi > 3 (n > 5). The points inside the caustic show multiply-imaged 
sources for which extra images are undetectable; the small (large) points de- 
note a flux ratio threshold / < 0.1 (f < 0.01). The inset shows a close-up of 
the tip of the inner caustic. 



iHolder & SchechteJ2003l) . and it is generally required for fit- 
ting obser ved galaxy-mass (M ~ 10 12 Mf^) strong lens sys- 
tems (e.g., ESonetaDIIlS IME^Maalllll). The lens 
potential associated with shear is 



1 

--fl 2 7Cos2(6>- 



6L), 



(7) 



where 7 is the dimensionless shear amplitude, 9~ ( is the shear 
direction. Shears of 7 ~ 0.05-0.1 are common for galaxy- 
mass lenses, and shears of 7 ~ 0.2-0.3 a re possible for lens 
galaxies lying in dense environments (e.g . ,|KeetOTi^talJl99J 



i Witt & Madfl997llKundic et alJll997albUFischer et alJU998 : 
kneib et alJ2000tlHolder & SchechteJ2003l) . 

3.2. Parameter Dependences for a Single Lens 

To begin to understand isothermal lenses, we show the 
source plane for a sample lens with ellipticity e = 0.5 in Fig- 
ure^ Singly-imaged sources with magnification /i > 3 occur 
only in a region just outside the caustics and near the minor 
axis of the radial caustic (which corresponds to the major axis 
of the galaxy density distribution). In this example with Ein- 
stein radius R e [ n =1", the image separations are larger than 
HST resolution, so the only systems that have only one de- 
tectable image are those whose extra images are too faint. 
Most of these are small flux ratio doubles, corresponding to 
sources that lie in two regions: just inside the outer radial 
caustic, where the secondary image is very faint; or near the 
inner tangential caustic, especially near the cusps along the 
major axis, where the primary image is highly magnified. In 
this example the two regions merge together when the flux 
ratio threshold for missing the second image is / < 0.1 (the 
small points in Figure^, but remain distinct when the thresh- 
old is / < 0.01 (the large points). 

We see that for isothermal galaxies with e = 0.5, the cross 
section for magnified systems with undetectable extra im- 
ages appears to be larger than the cross section for magnified 
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FIG. 2. — Maximum singly-imaged magnification for an isothermal sphere 
with shear (upper panel) and for an isothermal ellipsoid (lower panel). The 
maximum magnification becomes infinite at 7 = 1/3 or e = 0.606 (see the 
Appendix for details). 



singly-imaged systems. In the next two subsections we quan- 
tify this result carefully by computing the magnification cross 
sections for both cases. 

3.2.1. Singly-imaged magnification distributions 

For most isothermal lenses (unless the ellipticity or shear 
is large), there is a finite upper bound on the singly-imaged 
magnification. The bound is derived in the Appendix, and 
shown in Figure |2] While spherical models can never pro- 
duce a singly-imaged magnifica tion larger than Urnm = 2 (in 
the absence of microlensing; see iWvithe & Loe b 2002), non- 
spherical models can in principle produce much larger magni- 
fications. However, the effect is not likely to be very dramatic 
in practice: for a typical shear 7 ~ 0.1 or ellipticity e ~ 0.3, 
Mmax is still less than 3. 

The full magnification cross sections are shown in Figure[5] 
for various values of the parameters. Since all of the physi- 
cal parameters — the lens galaxy mass and redshift, and the 
source redshift — are contained in the Einstein radius, the 
dimensionless cross section A (pL)/R\ in depends only on ellip- 
ticity and shear so the parameter space we must study is small. 
Panels (a) and (b) show that shear and ellipticity increase not 
just fj, max but the whole high-magnification tail. However, the 
cross section for high magnifications is small; even in models 
with large shear or ellipticity, the area in the source plane with 
\i > 10 is less than 0.1% of the area with /i > 1.5. 

Shear and ellipticity are not mutually exclusive, so panel (c) 
shows what happens when we include both; we fix the ampli- 
tudes to typical values 7 = 0.1 and e = 0.3 and vary the angle 
between them. The effects are largest when the shear and 
ellipticity are aligned, because they combine to increase the 
quadrupole moment of the lens potential; and the effects are 
smallest when they are orthogonal because their quadrupoles 
partially cancel. When averaged over angle, we expect the 
combination of shear and ellipticity to produce a modest in- 
crease in the area with modest magnifications. 



3.2.2. Multiply-imaged magnification distributions 

Sample multiply-imaged magnification distributions are 
shown in Figure |4] There are various unusual features 
that can be understood with the help of the correspond- 
ing source plane shown in Figure [2 Without imposing 
any single-detectable-image criterion (SDIC), the magnifi- 
cation cross section shows a kink at the minimum magni- 
fication for quads, and the high-magnificat ion systems are 
dominated by quads (see. e.g..lBlandford & Kochanekl l 19871 
Schn eider etai1fT992l iFinch et alJl2002F SDICs remove the 
vast majority of quads, however, because the extra images in 
quads tend to be fairly bright. The only quad sources that sur- 
vive the cut lie extremely close to the caustic, with image con- 
figurations dominated by a very bright and very close pair of 
images. In this example, when the flux ratio SDIC is / < 0.1 
the magnification distribution is a smooth curve. However, 
when the flux ratio SDIC is / < 0.01 the distribution breaks 
up into two separate populations, with the low-magnification 
population lying just inside the radial caustic, while the high- 
magnification population lies just outside the cusps of the tan- 
gential caustic (see Figure^. (The image separation SDIC is 
unimportant here because the Einstein radius is larger than the 
HST resolution.) 

Having understood the general features, we can now exam- 
ine how the multiply-imaged magnification distribution de- 
pends on ellipticity, as shown in Figure [5] In the absence of 
SDICs, ellipticity raises the high-magnification tail; for ex- 
ample, the cross section for /i > 10 is increased by a factor 
of ~2 for e > 0.5. This case will apply to low-mass halos 
where the Einstein radius is small enough that all image sep- 
arations are unresolvable. When the flux ratio SDIC is im- 
portant (when the halo mass is large enough that the image 
separations would be resolved), ellipticity has a much more 
dramatic effect. In the spherical case there are no magnifica- 
tions jj, > 2/(1-/), where / is from the SDIC. Introducing 
ellipticity creates a population of high-magnification sources 
lying just outside the cusps of the tangential caustic. As e 
increases, that population grows and merges with the popu- 
lation of lower-magnification sources lying inside the radial 
caustic (as in Figure |4j}). Finally, as the ellipticity grows 
to e > 0.606 the cusp of the tangential caustic pierces the 
radial caustic, so the high-magnification region just outside 
the cusp becomes associated with singly-imaged rather than 
multiply-imaged systems. This explains why the multiply- 
imaged cross section curve changes shape, but it is not very 
important in practice because such large ellipticities are rare. 

The details clearly depend on the ellipticity and the SDIC, 
but the most important result is more general: the ordering 
of the curves in Figure [5] The dominant source of magnifi- 
cations /i > 5 should be unresolvable small-separation lenses 
produced by low-mass halos, followed by systems with mul- 
tiple images where the extra images are too faint to be de- 
tected. The contribution from true singly-imaged systems is 
generally not as important, except when the ellipticity is large 
(e > 0.6), and even then the cross section is quite small. 

Shear has similar effects on the magnification cross sec- 
tions, because like ellipticity it increases the quadrupole mo- 
ment of the lens potential and makes the tangential caustic 
larger. We do not show the cross sections for different shears, 
because the results appear very similar to those displayed in 
Figure|5] 

3.3. Averaging over Ellipticity and Shear 
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log /x log /x log ix 

FIG. 3. — Singly-imaged magnification distributions for isothermal models. The area A(p) where the magnification is greater than fi is expressed in units of 
^eiir ( a ) Effects of shear, (b) Effects of ellipticity. (c) Effects of ellipticity and shear together; the ellipticity and shear are fixed at e = 0.3 and 7 = 0.1, and we 
vary the angle between them. For reference, the dashed curve shows a model with e = 0.3 and no shear, while the dotted curve shows a model with 7 = 0. 1 and 
no ellipticity. 
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FIG. 4. — Multiply-imaged magnification distributions for an isothermal ellipsoid with ellipticity e — 0.5. The different line types denote different image 
configurations (doubles, quads, or all lenses). The different panels show different criteria for having only a single detectable image; panels (b) and (c) correspond 
to the detection limits in HST observations of the z ~ 6 SDSS quasars by Richards et al. 1 2004). 



In order to obtain the overall probability distribution for the 
magnification, we next average over realistic distributions of 
ellipticity and shear. For the distribution of shear amplitudes, 
we use the model derived by Holder & Schechter ( 20031) f° r 
the environments of early-type galaxies in A^-body and semi- 
analytic models of galaxy formation; they find a lognormal 
distribution with median 7 = 0.05 and dispersion <r 7 = 0.2 dex. 
We use random shear directions. For the ellipticity distribu- 
tion, we use data on the shapes of observed early-type galax- 
ies. 7 iTctrgensen. Franx. & Kiaergaardl Jl995) give ellipticities 
for 379 E and SO galaxies in 11 clusters, including Coma. 
The distribution is broad, with mean e = 0.31 and dispersion 
a e = 0.18. We average over more than 1000 random combi- 
nations of ellipticity and shear. 

Figure|6]shows the resulting cross sections. Singly-imaged 
systems are important only for low magnifications (/1 < 2). 
If the Einstein radius is small and the lenses are unresolved, 

7 The data give the shape of the light distribution, while what we need 
is the shape of the mass distribution. The mass and light shapes may not 
be correlated on a case-by-case basis, but for our purposes it is sufficient to 
assume that their distributions are similar (see Rusin & TegmarU200"II) . 



then multiply-imaged systems dominate at /x > 1.8; much of 
the relevant cross section comes from quads (as indicated by 
the difference between the solid and dashed curves). When 
the flux ratio SDIC applies, multiply-imaged systems domi- 
nate at (i > 2.5 for a SDIC of / < 0.1; for fi < 10 most of 
this cross section comes from doubles, while for higher mag- 
nifications there is a significant contribution from quads. For 
a flux ratio SDIC / < 0.01, multiply-imaged systems dom- 
inate at [i > 10, and most of the cross section comes from 
doubles. The overall conclusion is that for isothermal halos, 
most large magnifications /i > 10 will correspond to multiply- 
imaged systems where the extra image are not detectable (ei- 
ther unresolved or faint). 

We have explicitly examined ellipticity and shear, but 
early-type galaxies are also observed to have small oc- 
topole moments in their light distributions. We have re- 
peated our analysis using the ellipt icity and octopole d istri- 
butions from the galaxy samples of Bender et al. ( 1982) and 
Saglia riSender. & Dresslerl (119931) . and confirmed that our re- 
sults are not very sensitive to changes in the ellipticity distri- 
bution or to the addition of octopole terms. 
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FIG. 5. — Effects of ellipticity on multiply-imaged magnification distributions for isothermal models. In each panel, the three solid curves show different 
criteria for having a single detectable image; the top curve shows all systems, the middle curve shows A9 < 0"1 or / < 0.1, and the bottom curve shows 
A9 < 0"3 or / < 0.01. For comparison, the dotted curve shows the appropriate singly-imaged magnification distribution (see Figurel3*l. 




FIG. 6. — Magnification cross sections after averaging over ellipticity and 
shear. The dotted curve shows the singly-imaged cross section. The solid 
and dashed curves show the multiply-imaged cross sections. The top pair 
of curves includes all multiply-imaged systems (regardless of the number of 
detectable images). The middle pair corresponds to the following criteria for 
having a single detectable image: AS < 0. 1 i? C m or / < 0. 1 ; while the bottom 
pair corresponds the criteria A9 < 0.3i? e m or f < 0.01. 



4. NFWHALOS 

Another common and useful lens model is the Navarro- 
Frenk- White (NFW) profile, which describes halos produced 
in TV-body simulations. The NFW profile is thought to de- 
scribe systems that are dominated by dark matter at all radii: 
massive cluster halos, and perhaps low-mass dwarf halos as 
well. In this section we study the ability of NFW lenses to 
produce magnification without detectable multiple images. 



4.1. Definitions 
The NFW profile has the form 

Ps 



(r/r s )(l + r /r s ) 2 



(8) 



where r s is a scale radius and p s is a characteristic den- 
sity. There has been debate about whether the inner den- 
sity profile of simu lated clusters rea l ly asymptotes to the 
p oc r ~ l form (e.g., |NaAarro£Lalll^ 9^lF ukushige & Makino 
' iMooreet al.1 119991: LTing & Sutol 1200(1 iPoweretalJ 
lFukusMgeet^iTB004l) . and whether suc h a density 



1997; 
2003 



cusp is consistent with observed clusters (e.g., Tvsonetal. 
[I^ISmith et alJ200UlEttori et al.l2002UKelson et alJ2002t 
ISand et al.1 120021 Lewis et al.1 120031 ISand et al J 12004 . Our 
main need is for a model other than the isothermal ellipsoid 
that we can apply to massive clusters and low-mass dwarf ha- 
los. For this purpose the NFW model is standard and suffi- 
cient, and exploring a larger family of models (such as gener- 
alized NFW) is beyond the scope of this paper. Besides, with 
appropriate normalizations generalized NFW profiles lead to 
lens statistics that are no t so sensitive to the inner profile slope 
jKeeton & MadavJ2001l) . 

The NFW prof ile has projected surface mass density 
(Bartelma nnl 199611 



k(R) = 2k s 



\-F{R/r s ) 
(R/r s ) 2 -l 



(9) 



where k s = p s r s /S C nt is a dimensionless lensing "strength" 
parameter, and the function F(x) is: 



f (1 -jc 2 )-'/ 2 tanh"'(l -x 2 ) 1 / 2 x < 1 
F(x) =<l x =\ 
l(x 2 -l)- 1 / 2 tan- 1 (x 2 -l) 1 / 2 jc > 1 



(10) 



We obtain an elliptical NFW model by replacing R — + (x 2 + 
y 2 /q 2 y/ 2 in the surface mass density, where q < 1 is the pro- 
jected axis ratio, and the ellipticity is e = l-q. The lens- 
ing properties of an elliptical NFW model can be co mputed 
with a set of 1-d numerical integrals (Schramm 1990; Keetonl 
I2001bh . 8 

8 It is possible to obtain an analytic NFW model by putting the ellipti- 
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FIG. 7. — Sample lensing strength k, (upper panel) and halo scale radius 
r s (lower panel) versus halo mass, for a lens at redshift zi - 1 and source 
at redshift Zs = 6. The solid curves show results for halos with the median 
concentration, using the median c(M) relation from Bullock et al. 1 2001); the 
dotted curves indicate the 1<t range due to the scatter in the M-c correlation, 
a c = 0.14 dex at fixed mass. 



NFW profiles appear to form a two-parameter family spec- 
ified by p s and r s , or equivalently by the virial mass M and a 
concentration parameter c = r v i, /r s . (The virial radius r v i r can 
be given as an explicit function of M and c.) In fact, the two 
parameters are correlated, and different mode ls for m edian 
relation and scatter h ave been proposed (e.g.. [N avarro et al.l 
f!997llEke et alJ200UlBuTlock et alJ200!llJing & Sutol2002l) . 
For our purposes, the important result is that the scale length 
r s and lensing strength k s are both correlated with the halo 
mass, as shown in Figure In this section we can express 
cross sections in units of r 2 so we need not examine the r s 
dependence explicitly; but we do need to examine the depen- 
dence on k s . 

4.2. Parameter Dependences for a Single Lens 

To begin to understand NFW lenses, in Figure [8] we show 
the source plane for a sample lens with strength k s = 0.168 
(corresponding to a median 10 14 M Q halo at zi = 1, see Fig- 
ure and ellipticity e = 0.1. The caustics are small, and 
the region of multiply-imaged systems that satisfy reason- 
able single-detectable-image criteria (SDIC) is smaller still. 
The reason is that the image separations for massive halos are 
large, so the important SDICs involve the flux ratios. NFW 
halos generally produce large magnifications, and the range 
of magnification ratios is not very broad. As a result, the 
only sources that survive the SDICs lie very near the caustics, 
where one of the images has a very large magnification. 

By contrast, there is quite a large region of the source plane 
where sources are singly-imaged but have large magnifica- 
tions. Thus, for NFW halos it appears that singly-imaged sys- 

cal symmetry in the potential r ather than the density (e.g., Golse & Kneihl 
2002; Meneahetti, Bartelmann, & Moscardini 2003). However, /V-body sim- 
ulations suggest that it is the density rather than the potential that has ellip- 
soidal symmetry (or more generally, triaxiality; e.g., Jing & Suto 2002). We 
find that working with an elliptical density and using numerical integrals is 
not a major hindrance. 



tems will be more important than multiply-imaged systems 
for producing large magnifications with a single detectable 
image. In the next two subsections we quantify this result 
carefully. 

4.2.1. Singly-imaged magnification distributions 

Even spherical NFW lenses are complex systems where the 
lens equation is transcendental, so they must be studied nu- 
merically. Figure [5] shows the maximum singly-imaged mag- 
nification as a function of the strength k s . Spherical NFW ha- 
los can apparently produce large magnifications without mul- 
tiple imaging, especially when the lensing strength (or halo 
mass) is small. It has been known that for multiple imaging at 
fixed splitting angle NFW lenses tend to produce smaller cross 
sections but larger magnifications than i sothermal lenses (see 
iBlandford & Kocha nekTT987l iKnudson et al!2001l) . Now we 
see that the association between NFW lenses and high magni- 
fications extends to single imaging as well. 

The full magnification distributions are shown in Fig- 
ure [H)] Panel (a) shows that for spherical halos, increasing 
the strength k s decreases /i max but increases the overall cross 
section. Note that the figure shows the cross section in units of 
r 2 ; expressing the area in physical units like steradians would 
increase the distance between the curves, since the strength is 
correlated with r s (see Figure^}. 

Panels (b) and (c) in Figure^|show that shear and elliptic- 
ity do not dramatically affect the magnification distributions 
for NFW halos. This stands in contrast to the case for isother- 
mal halos. The difference is that the magnification distribu- 
tion for spherical NFW halos already extends to high magnifi- 
cations, so any increase due to ellipticity or shear is relatively 
more moderate. The departure from spherical symmetry does 
raise the tail to very high magnifications, but that effect is 
not very sensitive to the degree of asymmetry (shear or ellip- 
ticity). Interestingly, ellipticity appears to lower slightly the 
cross section for moderate magnifications. Finally, FigurelTUtl 
shows that allowing a combination of ellipticity and shear has 
little effect other than extending the high-^i tail of the distri- 
bution. While this is shown explicitly only for fixed e = 0.3 
and 7 = 0.1, we expect this generic feature to hold for other 
combinations. 

4.2.2. Multiply-imaged magnification distributions 

Figure ^2 shows sample multiply-imaged magnification 
distributions for spherical NFW halos with different values 
of the lensing strength k s . Only the case with no regard for 
the number of detectable images is shown, because as we saw 
in Figure [8] the cross section for NFW lenses that satisfy the 
flux SDICs is extremely small. When the strength is low, the 
cumulative cross section A(p) is flat out to p > 100, indicat- 
ing that all magnifications are larger than 100; however, the 
cross section is very small. As k s increases, the point at which 
A(p) begins to decline moves to the left, indicating that the 
minimum magnification decreases; and the cross section in- 
creases. For all but the most massive and concentrated halos, 
the multiply-imaged magnification cross section is consider- 
ably smaller than the singly-imaged cross section, except at 
the very highest (and rarest) magnifications. 

A better way to compare the singly- and multiply-imaged 
cases is to take the ratio of the cross sections, as shown in Fig- 
ure [21 Decreasing k s increases the ratio dramatically, so for 
moderate- to low-mass NFW halos the singly-imaged magni- 
fication cross section can be orders of magnitude larger than 
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arcseconds arcseconds 

FIG. 8. — Source plane for an NFW halo with strength k s = 0.168, scale radius r s = 167 /T 1 kpc = 29"8, and ellipticity e = 0.1. The curves show the 
caustics. The small (large) points outside the caustics indicate singly-imaged sources with magnification /x > 3 (/i > 10). The small points inside the caustic 
show multiply-imaged sources for which extra images are undetectable; the criteria for having only a single detectable image are an image separation A8 < O'.'l 
or flux ratio / < 0. 1 for the small points, and AS < 0"3 or / < 0.01 for the large points. The left panel shows a large region of the source plane, while the right 
panel shows a close-up of the multiply-imaged region. 
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FIG. 9. — Maximum singly-imaged magnification versus lensing strength 
for spherical NFW lenses. 



the multiply-imaged cross section. This quantifies the state- 
ment that for magnification by NFW halos, singly-imaged 
systems are vastly more important than multiply-imaged sys- 
tems. 

Figure ^] shows that ellipticity has little effect on the 
multiply-imaged magnification distribution. Large shears can 
in principle increase the cross section. However, it is not 
clear that massive NFW halos could experience such large 
shears. Shears of 7 ~ 0.2-0.3 typically occur when the lens 
is a galaxy within a clusters, where the cluster serves as the 
environment that produces the shear. The filamentary struc- 
ture typical around clusters represents a very different envi- 
ronment, for which the shear distribution is not well known. 
Fortunately, this uncertainty is not important for our results 
because shear and ellipticity have such modest effects. For 
simplicity, in the rest of the paper we consider only spherical 
halos when computing multiply-imaged magnification distri- 



butions for NFW halos. 

5. A REALISTIC HALO POPULATION 

Having understood two fiducial halo models, we now com- 
bine them into a realistic population of galaxies and clusters. 
We compute the overall probability for magnification with ad- 
ditional detectable images, and use it to evaluate the hypoth- 
esis that the four z ~ 6 SDSS quasa rs are highly amplified by 
lensing. After defi ning the model (§ 15. 1> . we present our gen- 
er al re sults (§ I5.2> and then apply them to the SDSS quasars 
(§ 15. 3> . We end with a discussion of some systematic effects 
in our analysis (§ 15. 4> . 

5.1. The Model 

The total probability distribution for fi comes from integrat- 
ing the cross section over an appropriate halo population, 

P ^z s )= fdvj dM^A(fx;z s ,zi,M)- (11) 

The first integral is over the comoving volume between the 
observer and source. The second integral is over the co- 
moving halo mass fu nction dn / dM\ we ado pt the theoretical 
mass function from Sheth & Tormen ( 1999). Finally, A(p) is 
the cross section computed above (expressed in steradians), 
which depends on the source redshift z s , the lens redshift zi, 
the lens halo mass M, and also on the lens model (isother- 
mal or NFW). By using the appropriate cross section, we 
can compute the probability for singly-imaged or multiply- 
imaged magnifications by any of the halo populations we have 
considered. 

Implicit in eqn. il Q is an average over appropriate ellip- 
ticity and shear distributions. For isothermal halos we use 
t he re sults after averaging over ellipticity and shear, from 
§ 13.31 Since the results for NFW halos are not very sensi- 
tive to ellipticity and shear, we simply use a fixed ellipticity 
e = 0.1 for the singly-imaged case (in order to pick up the 
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FIG. 10. — Singly-imaged magnification distributions for NFW models. The area A(fi) where the magnification is greater than fi is expressed in units of 
(not ^f?j_). (a) Effects of the lensing strength k s ; the curves show re s = 0.08, 0.10, 0. 14, 0.20, 0.30 from bottom to top. (b) Effects of shear; the curves show 
7 = 0,0.1,0.2,0.3,0.4 from bottom to top. (c) Effects of ellipticity; the curves show e = 0,0.2,0.4,0.6,0.8 from top to bottom, (d) Effects of ellipticity and shear 
together; the curves show A6 = 0,30,60,90 deg. The ellipticity and shear are fixed at e = 0.3 and 7 = 0.1. For reference, the dotted curve shows a spherical 
model. In panels (b)-(d) the lensing strength is fixed at K s = 0. 168 (the median value for a 10 14 h Mq halo in Figure 171. 
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FIG. 1 1 . — Magnification distributions for spherical NFW halos with different values of the lensing strength k s . The solid curves show multiply-imaged mag- 
nification distributions, r egar dless of the number of detectable images. For comparison, the dotted curves show the corresponding singly-imaged magnification 
distributions from Figure Pol 




cation cross sections for NFW lenses. The curves correspond to k„ = 
0.08,0.10,0.14,0.20,0.30 from top to bottom. 



high-magnification tail; see Figure [Tot, and we use spherical 
models for the multiply-imaged case. 

We consider a model with at least two different halo popula- 
tions. The most massive halos, corresponding to clusters and 
groups of galaxies, are treated as NFW halos; while halos cor- 
responding to galaxies are modeled with isothermal pr ofiles. 
In this increasingly standard model llFlores & PrimacklH99^ 



Keeton 1998; Porciani & Madau 2000; Kochanek & White 
120011 Li & Ostriker 2002J), the difference between clusters 
and galaxies is usually attributed to baryonic cooling: in mas- 
sive halos the baryons have not had time to cool so the systems 
retain their initial NFW form; while in lower-mass halos the 
gas has cooled and condensed into the center of the system, 
and created a mo re concen trated tot al (baryons + dark matter) 
mass profile (e.g.. lBlumenthal et alJl986llKochanek & White! 
2001). Thus, the transition between clusters and galaxies 
is characterized by a mass scale M c i us such that halos with 
M < M c ius (M > M c i us ) have a cooling time shorter (longer) 
than the age of the universe and correspond to galaxies (clus- 
ters). 

There may be a third population as well, namely low- 
mass "dwarf" halos with NFW profiles. Theoretical halo 
mass functions rise much more steeply than observed galaxy 
luminosity functions, leading to the speculation that there 
may be a substantial population of underluminous low- 
mass halos. Mechanisms such as feedback or reionization 
might have suppressed bary onic cooling and sta r formation 
in lo w-mass systems (e.g., IP ekel & Silkl 119861 lEfstath ioul 
1 19921 iNavarro & Steinm etzlll997t iThoul &" Weinbergf Tl99a 
Bullock et al. 2000; Springel & Hernquist 2003), leading to 
halos that are dark and retain their initial NFW form. The ap- 
parent dearth of small-separation lens systems implies a tran- 
sition from isothermal galaxies back to NFW dwarfs around 
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FIG. 13. — Magnification distributions for NFW halos with different ellipticities (upper panel) and shears (lower panel), with no regard for the number of 
detectable i mag es. The lensing strength is fixed at k s = 0.168. For comparison, the dotted curves show the appropriate singly-imaged magnification distributions 
from Figure fTol 
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FIG. 14. — Net singly-imaged magnification probability distributions for sources at redshift Zs = 6- P(n) is the cumulative probability of having a magnification 
greater than fi. The solid curves show the total probability, while the short-dashed, dotted, and long-dashed curves show the contributions from dwarf halos, 
normal galaxies, and clusters, respectively. The different panels correspond to different assumptions about the dwarf/galaxy transition (see text). 



a mass of Md wa rf ~ 10 12 M Q iLi & Ostrikerj|2003t lMall2003t 
IKuhlen et alJ20 04). 

Following iMal fe003). we summarize the model by intro- 
ducing a function /sis(M) that describes the fraction of ha- 
los of mass M that have isothermal profiles (and the rest are 
NFW). We use the function: 



/sis(M) = 



exp 



(logM-logM dus r 



M < M dwarf 

Md„ ar f < M < M dus 

M > M dus 



112) 



To det ermine the parameter values, we follow IKuhlen et alJ 
(2004) and fit the model to the observed image separation 
distribution from the Cosmic Lens All -Sky Survey (CLASS; 
iMvers et alJl2003l iBrowne et alJl2003l) . 9 The formal best-fit 

9 Our quantitative results differ slightly from the fiducial results of 



model has a c \ us =0.13 and logM c i us = 13.22. However, there 
is a degeneracy between cr c i us and M c i us such that a model 
with a sharp transition (cr c i us = 0) and logM c i us = 13.38 fits 
almost as well. We have verified that the two models give 
indistinguishable results, so we report results only for the 
model with a sharp transition. The cluster/galaxy transition 
mass inferred from lensing ag rees well with estimates based 
on cooling arguments ("e.g.. iKochanek & White! 1200 lb . As 
for the galaxy/dwarf transition, current data lack any ability 
to determine whether it is smooth or sharp, so for simplic- 
ity we use only a sharp transition. The formal best-fit model 
has logMd„aif = 12.5, although in fact onl y an upper limit on 
Mdwarf is reliable (see Kuhlen et al. 2004). This mass is too 



IKuhlen et all 1 2004 ) because we now use the 


Sheth & Tormen 1 1999) 


function rather than that from Jenkins et alJ 


2001), and we use cr% = 



0.90 



rather than cr 8 = 0.74. 
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high to be explained by reionization feedback, so it might in- 
dicate either some other kind of feedback or some peculiari- 
ties in the data (perhaps unknown incompletenesses at small 
separations, or simply small-number statistics). We therefore 
consider several different values for M dW arf- We take the tran- 
sition masses M c i us and Mdwarf to be independent of redshift. 
This assumption may seem objectionable, but we show below 
that the specific value of M c ] us has little effect on our results, 
and we explicitly consider systematic uncertainties associated 
with M dwarf . 

5.2. Main Results 

Figure[2]shows the net singly-imaged magnification distri- 
bution, for our fiducial model. For sources at redshift z s = 6, 
magnifications of \i ~ 1 .6 occur at the percent level, and the 
distribution drops quickly. The probabilities for /i > (2,5, 10) 
are (0.3%, 9 x 1 -5 ,1.1 x 10" 5 ). In other words, significant 
singly-imaged magnifications are rare. 

There are nevertheless several qualitative results that are in- 
teresting and instructive. First, for p, > 1.8 the singly-imaged 
magnification probability is dominat ed by clusters. Th i s resul t 
seems at first glance to contradict Comerford et al. (2002), 
who found that massive NFW halos had negligible impact 
on the probability for lensing magnificatio n, but the appar- 
ent discrepancy is easily explained. Comerf ord et alJ consid- 
ered lensing with undetected multiple imaging, while we are 
considering the complementary case where there is only a sin- 
gle image. NFW halos are much more efficient at producing 
highly magnified single images than multiple images (see Fig- 
ure I12> . which is why we find a much stronger effect. The 
important implication of this result is that if there is a sig- 
nificant singly-imaged magnification (/i > 5) then the lensing 
object is most likely a cluster, and ought to be relatively easy 
to detect. Conversely, lensing with small to moderate singly- 
imaged magnifications (/i < 2 ) is dominated by ga laxies. This 
is consistent with the claim by Shiova et al. ( 2002) that one of 
the z«6 quasars is magnified by a factor /i w 2 by a fore- 
ground galaxy, with no obvious sign of a cluster. 

A second interesting feature of the singly-imaged magni- 
fication probability relates to the possible galaxy/dwarf tran- 
sition at low masses. As discussed above, current lens data 
hint at the transition but do not constrain it well, so we con- 
sider three possible cases: (a) a model with no transition; 
(b) a model with the transition at logMi„ ar f = 12.0; and (c) 
a model with the transition at logMd war f = 12.5, which per- 
haps seems high but is formally the best fit to current data 
(although it is not significantly better than the other models). 
Models (a) and (c) are extremes that bound the range of rea- 
sonable possibilities, and model (b) is a sample intermedi- 
ate model. Not surprisingly, Figure ^] shows that moving 
the transition changes the relative contributions of galaxies 
and dwarf halos to the probability; if the transition occurs at 
logMdwarf 12.0, dwarfs can conceivably contribute more of 
the probability than normal galaxies. But the important result 
is that changing Mi„ ar f affects the net probability by <10%, 
since the singly-imaged case is dominated by clusters. 

Figure[2]shows the multiply-imaged magnification proba- 
bility, for different assumptions about the dwarf/galaxy tran- 
sition and different single-detectable-image criteria (SDIC). 
In general, nearly all of the probability comes from galaxies. 
Clusters contribute < 10~ 5 of the probability if we consider all 
image configurations, and < 10~ 8 if we consider those with 
only a single detectable image; so we confirm the result from 
IComerford et alJ (12002) that clusters are essentially negligi- 
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NOTE. — Fractio n of lens systems m agnified by fi > 10 that lack extra im- 
ages detectable by Richards et al. 1 2004). Columns 2-3 give the fractions that 
are singly-imaged and multiply-imaged, respectively, and Column 4 gives 
the total. Here we give results for a model with the dwarf/galaxy transition at 
logMd war f = 1 1.0, in addition to the three models discussed in the text. 



ble for the multiply-imaged magnification probability. NFW 
dwarfs halos, if present in the model, likewise have a neg- 
ligible contribution to the probability (no more than ~ 10~ 7 
even in the model with logMd„ a rf = 12.5). The main effect 
of having dwarfs in the model is to reduce the number of ha- 
los that are isothermal galaxies, and hence to reduce the net 
multiply-imaged magnification probability. The reduction can 
be substantial when the SDICs are important. 

Perhaps the most interesting result from Figure ^] is the 
effects of the SDICs. Consider the probability of a magni- 
fication fi > 10, and suppose we have no knowledge of the 
presence or absence of additional images. Then we must use 
the multiply-imaged probability with no SDICs, and add the 
singly-imaged probability, which yields total probabilities of 
(2.1, 1.4, 1.0) x 10" 4 f or models (a), (b), and (c) respectively. 
However, suppose like Richards et al. (2004) we can rule out 
the presence of extra images down to either Ad = Of 1 and 
/ = 0.1, or A9 = 0."3 and / = 0.01. We should then use 
the lowest of the multiply-imaged curves in Figure ^] and 
again add the singly-imaged probability. A useful way to 
quantify the results is to give the fraction of systems mag- 
nified by u > 1 tha t do not have extra images detectable by 
Richards et al. (2004). This fraction is given by the ratio of 
the probability with SDICs to the probability without, or dif- 
ferences between the curves in Figure^J the results are given 
in Table For example, in the no-dwarf model, about 24% 
of highly magnified systems are multiply-imaged such that 
the extra images are undetectable, and another 5% are singly- 
imaged. In models with dwarfs, the fraction that are multiply- 
imaged without detectable extra images is much lower, so in 
total only 9-15% of highly magnified systems lack detectable 
extra images, and many of those are true singly-imaged sys- 
tems lensed by clusters. 

5.3. Implications for the SDSS quasars 

So far we have only considered a priori lensing probabili- 
ties, i.e., the bare optical depth for producing certain magni- 
fications. However, to compute the probability of finding a 
certain magnification in a real, flux-limited survey we would 
also have to fold in magnification bias to obtain a posteriori 
probabilities. Because the SDSS can probe only the steep, 
bright end of the quasar luminosity function at z ~ 6, magni- 
fication bias may be quite strong, but it is is extremely sen- 
sitive to the poorly-known LF shape. (In fact, the problem 
can be turned around so that lensing [or lack thereof] in the 
z ~ 6 quasar sa mple yiel ds constrains on the L F slope; see 
IComerford et al.ll2002l and lRichards et al.ll2004h Rather than 
making detailed but highly LF-dependent predictions of the 
a posteriori probabilities, we turn attention to probability ra- 
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FIG. 15. — Net multiply-imaged magnification probability distributions for sources at redshift z s = 6. In each panel, the upper solid curve shows the overall 
multiply-imaged probability (regardless of the number of detectable images), and the two other solid curves show the probability with different criteria for 
detecting only a single image (as indicated). The dashed curves show the contribution to each probability from galaxies alone (which is often indistinguishable 
from the total). The different panels correspond to diffe rent assumptions about the dwarf/galaxy transition. For comparison, the dotted curve shows the net 
singly-imaged magnification distribution from Figure fT4l 
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FIG. 16. — Similar to Figure [l5l but for a model in which all halos are spherical and there is no shear. 



tios and consider the following question: What is the ratio of 
the probability for /j, > 10 with observational SDICs, to the 
probability for fj, > 10 without SDICs? This is equivalent to 
the question: In a toy model with magnification bias so strong 
that all z~6 quasars are magnified by a factor /i > 10, what 
is the probability that one o f the quasars would h ave no addi- 
tional images detectable by Richards et al. (2004? 

In general, probability ratios should be the same for a pri- 
ori and a posteriori probabilities, so Tabled is exactly what 
we need. To summarize, if we are extremely optimistic about 
having strong magnification bias and about all low-mass ha- 
los having steep isothermal profiles, we can imagine that the 
probability of a z ~ 6 quasar being magnified by /i > 10 
without having additional images detectable by HST might 
be as high as 29%. If this is the case, we cannot rule out 
the possibility t hat on e of the four z ~ 6 quasars observed by 
iRichards et alJ fc004l) might actually be significantly magni- 
fied. However, even then the probability that all four are am- 
plified would be P = (0.29) 4 = 0.007. The actual probability 
is almost certainly much lower, because magnification bias is 
probably not as strong as in the toy example, and because as- 



suming that all low-mass halos have isothermal profiles prob- 
ably overestimates the lensing optical depth. In other words, 
we can rule out at more than 99.3% confidence the hypothesis 
that all four z ~ 6 quasars are amplified by more than a fac- 
tor of 10 — provided we can equate a priori and a posteriori 
probability ratios. 

The only possible problem with equating the probability ra- 
tios is if the SDSS is somehow biased against high-redshift 
lensed quasars. Suppose, for the sake of argument, that 
there were a total of 4/0.29 = 14 lensed and highly magni- 
fied z « 6 quasars in the SDSS, 10 of which have not been 
identified. Then having f our quasars tha t are m agnified but 
lack images detectable by Richard s' et alJ (|2004) would actu- 
ally be consistent with the no-dwarf model. In our models 
with dwarfs, the total number of lensed quasars needed would 
be ^25^45 (see Table 0. The key question is whether so 
many "missing" lensed high-redshift quasars could exist. The 
SDSS z > 5.8 quasar sample is selected on the basis of col- 
ors alone; the sample is basically /-dropouts, with some addi- 
tional color criteria to reduce contaminants ( Fan et alJl2000l 
1200 It 120031) . There is no requirement that the objects ap- 
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pear point-like in SDSS images, and hence no bias against 
multiply-imaged systems. The only remaining problem is if 
light from a lens galaxy could change the composite colors 
of a lensed quasar system enough that the system would not 
be selected as an /-dropout. Lens galaxies associated with 
z ~ 6 lensed quasars would be expected to l ie at redshifts 
1 < zi < 2 (see Fig. 1 of Comerford et al. 2002), and so would 
pro bably be too faint to s ignificantly change the colors (also 
see Wvithe & Loebl 120021) . While it would be interesting to 
quantify this effect more carefully, 10 it seems unlikely that 
the SDSS high-redshift quasar sample is highly biased against 
multiply-imaged systems. 

5.4. Corollaries and Systematics 

We have obtained our main result, namely answering the 
question of whether it is likely that the z ~ 6 quasars are highly 
magnified, but there are several corollaries worth mentioning. 
First, so far in this section we have allowed the galaxies to 
have ellipticity and shear. (We have made certain assumptions 
about the ellipticity and shear f or dw arf and cluster halos, but 
they are not so important; see § 15.11 ) It is interesting to repeat 
the analysis with all halos assumed to be spherical, to see how 
much ellipticity and shear affect the results. Figure [161 shows 
the results. In the no-dwarf model, neglecting ellipticity and 
shear reduces the probability of multiple imaging with fi > 10 
by 30-50%. The same holds in models with dwarfs for the 
total multiple imaging probability (with no selection effects). 
However, in models with dwarfs where the SDICs apply, the 
reduction is at least a factor of 5 and often much larger. This 
result can be understood as follows. Most of the multiply- 
imaged magnification probability comes from the isothermal 
halos called "galaxies" in our model. If there are no dwarfs, 
then most of the probability comes from low-mass galaxies 
that produce image separations too small to be resolved; so 
the flux ratio SDIC is unimportant, and the probability is not 
dramatically sensitive to ellipticity (see the upper curves in 
Figure |3- By contrast, if low-mass halos are dwarfs (which 
have negligible cross sections), then much of the probabil- 
ity must come from galaxies that produce image separations 
larger than the resolution. In this case the flux ratio SDIC 
plays an important role; and in Figure|5]we saw that this dra- 
matically reduces the cross section when the ellipticity is zero. 
The bottom line is that ellipticity and shear are a factor of ^2 
effect in no-dwarf models, but can be an order of magnitude 
effect in models with dwarfs. 

The second point is that there are several additional sys- 
tematic uncertainties that might be relevant. We have con- 
sidered three effe cts that have the m ost impact on strong 
lens statistics (see Kuhlen et alJ|2QQ 4j): fi) changing the halo 
m ass function f r om th at of Sheth & Tormen ( 1999) to that 
of iJenkins et alJ d20Qll) : (if) changing the scatter in the M- 
c correlation for NFW halos from 0.14 dex to 0.07 dex or 
0.21 dex; and (iii) varying the location of the cluster/galaxy 
transition by A(logM c i us ) = ±0.25. We find that changing the 
mass function has the strongest effect, and even that is only 
a ~10% change in the net magnification probability. In other 
words, our results appear to be robust to effects other than the 
question of whether low-mass halos are isothermal galaxies 
or NFW dwarfs. 

Finally, our discussion has been geared toward sources at 

10 An analysis that begins with the empirical correlation between im- 
age sep aration and lens galaxy luminosity (see the Appendix of R usin et af] 
2003a) would be reasonably straightforward. 
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FIG. 17. — Total magnification probability distributions, including both 
the singly-imaged and multiply-imaged contributions, as a function of source 
redshift is = 6,5,4,3,2,1 from top to bottom. Here we use the no-dwarf 
model to obtain the maximum possible probability. 

Z « 6, but for completeness in Figure El we show the net a 
priori magnification probability as a function of source red- 
shift. We use the no-dwarf model as a way to obtain an upper 
bound on the probability; for the multiply-imaged case we use 
the criteria A9 < 0."3 or / < 0.01 for detecting only a single 
image, although from Figurefl5h this choice is not so impor- 
tant. Reducing the redshift naturally reduces the probability, 
especially for z s < 3. It affects the whole distribution in the 
same way, so our general conclusions appear not to be highly 
sensitive to the source redshift. Pushing quasar and galaxy 
samples beyond z ~ 6 will not significantly increase the prob- 
ability for large lensing magnifications. 

6. CONCLUSIONS 

The problem of lensing magnification without multiple de- 
tectable images has a rich phenomenology. First, there is the 
case of true singly-imaged systems. Isothermal halos are not 
very efficient at producing highly magnified single-image sys- 
tems, but NFW halos are. Consequently, high singly-imaged 
magnifications are possible in principle, and they are mainly 
associated with massive (> 10 13 5 M Q ) halos corresponding to 
clusters of galaxies. An important implication for observa- 
tions is that, if there is no evidence for a cluster along the 
line of sight to a distant quasar, then it is unlikely that there is 
strong singly-imaged magnification. 

The second case is when there are multiple images but the 
extra images are not detectable, either because the image split- 
ting is too small to be resolved, or because the extra images 
are too faint. NFW halos are inefficient at producing multiple 
images, and when they do they rarely produce extreme flux ra- 
tios; therefore, clusters contribute negligibly to the multiply- 
imaged magnification probability. Instead, this case is domi- 
nated by galaxies and lower-mass systems (< 1O 13 M ). The 
probability is very sensitive to the inner density profile of 
these halos. If all low-mass halos have steep isothermal pro- 
files, then the probability is dominated by lens systems with 
image separations too small to be resolved by HST However, 
if low-mass dwarf halos have NFW profiles (and hence small 
cross sections), then the overall probability is dominated by 
lens systems where the extra images are faint. 

Our central quantitative result is that 9-29% of all lens sys- 
tems with magnifications /i > 10 lack additional detectable 
images. In a toy model where magnification bias is so strong 
that most or all z ~ 6 quasars are lensed, then we cannot rule 
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out the hypothesis th at one of the four quasars observed by 
iRichards et alJ fc004l) is magnified despite lacking extra im- 
ages. However, even in such an extreme model, the prob- 
ability that all four are magnified by a factor of 10 would 
still be no more than 0.7%, and is probably much lower. The 
only way to evade this argument is if the SDSS high-redshift 
quasar sample is somehow biased against quasar lens systems 
with multiple bright images. That seems unlikely, although 
a detailed analysis of whether light from a lens galaxy could 
cause the /-dropout selection technique to miss a lens system 
is needed to answer this question definitively. 

Incide ntally, we can comme nt on the two different criteria 
used by Richar ds et alJ ll2004l) to search for companion im- 
ages to the z ~ 6 quasars. Richards et al. were able to rule out 
extra images down to a flux ratio / = 0.01 at image splittings 
greater than 0."3, or down to a less stringent flux ratio / = 0. 1 
for smaller image splittings down to 0"1. We find that for 
the no-dwarf model, the two criteria give similar probability 
results. By contrast, for models with dwarfs, the criteria with 
more stringent flux bounds provide stronger constraints on the 
lensing probabilities. In other words, for the purpose of de- 
termining whether distant quasars are magnified by lensing, it 
is more valuable to aim for more dynamic range than to push 



for resolution much better than ^0"3. 

To summarize, if the SDSS high-redshift quasar sample 
is not highly biased against multiply-imaged quasars, then 
it is quite im probable that all four quasars observed by 
IRichards et alJ J2004I) are highly magnified. In that case, the 
quasars can be taken as good evidence for the presence of 
billion-M black holes in the young universe. Explaining 
such black holes is a challenge for black hole growth mod- 
els, whose solu tion may involve a need for super-Eddington 
accretion (Haiman 2004). In other words, in the case of high- 
redshift quasars, a lensing "null result" actually makes the ob- 
jects even more interesting. 
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APPENDIX 

MAXIMUM SINGLY-IMAGED MAGNIFICATION FOR ISOTHERMAL LENSES 

With isothermal lenses there are three simple cases in which the maximum singly-imaged magnification can be obtained 
analytically; the bounds were are shown in Figure|2]and are derived here. First, for a simple isothermal sphere the magnification 
of a singly-imaged source (u > 7? e ; n ) is \i = 1 +R e i n /u, so the maximum magnification for singly-imaged sources is /i max = 2. In 
fact, in this case the full singly-imaged magnification cross section can be derived analytically, 
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Next consider an isothermal sphere with external shear. The magnification as a function of position is 
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(We are now working in a coordinate system aligned with the shear, so 9 1 = in these coordinates.) The radial caustic is a circle 
with radius R e i„. It maps to a curve in the image plane called the 1-2 transition locus, which marks t he transition from single 
images to images that are part of a two-image system, and which can be written in polar coordinates as (Finch et al. 2002) 
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For 7 < 1/3 the radial caustic completely encloses the tangential caustic, so all sources outside the radial caustic and all images 
outside the 1-2 transition locus are singly-imaged. In this case, the maximum singly-imaged magnification occurs at 6 = on the 
1-2 transition locus, 

Mmax = n _ 2 V - . • (7 < 1 /3) (A4) 

(1 -37)(l+7) 

For 7 > 1/3 the tangential caustic pierces the radial caustic to form a naked cusp (e.g.. lSchneider et all 1 9921 iFinch et all2 002). 
Sources just outside a naked cusp are singly-imaged but can have arbitrarily large magnifications, so /i max — > 00 in this case. 
Finally consider an isothermal ellipsoid. The magnification as a function of position is 
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We cannot compute the full 1-2 transition locus, but we can compute where it intersects the x and y axes: 
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The critical value of the axis ratio at which the tangential caustic pierces the radial caustic to form a naked cusp is found by 
solving 



2 tan 
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= o, 



(A8) 



whose solution is q = 0.394, or ellipticity e = 0.606. For q < 0.394, the presence of a naked cusp again causes /i max — > oo. For 
q > 0.394, the maximum singly-imaged magnification again occurs at 9 = on the 1-2 transition locus and has the value 
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